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00 . Measurements of flux-tubes generated by sources in different representations 



of SU(2) and U(l) lattice gauge theory in three dimensions are reported. Heavy 



' "quarks" are considered in three representations of SU(2): fundamental (j = 1/2), 

O ■ adjoint (j = 1), and quartet {j = 3/2). Wilson loops are used to introduce a 

. static quark-antiquark (QjQj) pair. Several attributes of the fields generated by 

the QjQj pair are measured. In particular, the first direct lattice measurements of 
the flux-tube cross-section Aj as a function of representation are made. It is found 
p ' that Aj ~ constant, to within about 10% (a rough estimate of the overall quality 

of the data). The results are consistent with a connection between the string ten- 
sion aj and cross-section suggested by a simplified model of flux-tube formation, 
aj = g^j{j + l)/{2Aj) [where g is the gauge coupling], given that the string tension 
^ I scales like the Casimir j{j + 1), as observed in previous lattice studies in both three 

' and four dimensions (and confirmed here up to the quartet representation). These 

results can be used discriminate among phenomenological models of the physics un- 
derlying confinement. Flux-tube measurements are also made in compact QED3, 
which exhibits electric confinement due to magnetic monopole condensation. Singly- 
and doubly-charged Wilson loops are considered. The string tension is found to scale 
like the squared-charge, and the fiux-tube cross-section is found is be independent of 
the charge, to a good approximation. The results of these three-dimensional SU(2) 
and U(l) simulations taken together lend some support, albeit indirectly, to a con- 
jecture that the dual superconductor mechanism underlies confinement in compact 
gauge theories in both three and four dimensions. 
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I. INTRODUCTION 



Flux-tube formation provides an attractive description of confinement in quantum chro- 
modynamics (QCD). In a simplified flux-tube picture of a very heavy quark-antiquark pair, 
color-electric field lines running between the quarks are assumed to be "squeezed" into a 
cylinder whose cross-section is independent of the quark separation R (ignoring end effects), 
resulting in linear confinement. 

Consider heavy "quarks" in an arbitrary representation of the gauge group, which we 
take to be SU(2) for convenience. The color-electric field Ej between quarks in the j-th 
representation is determined, in an Abelian approximation, by Gauss' Law [|I| EjAj = gQj, 
where the quark color "charge" Qj is related to the group Casimir, Q"^ = j{j + 1), g is the 
gauge coupling, and Aj is the cross-section of the flux tube. In this simplified model, Ej is 
assumed to be constant across the flux-tube cross-section. In the case of non-fundamental 
representation sources, R is also assumed to be below the threshold for fission of the tube. 
The interaction energy of the system Vj^^{R) = ^EjAjR = ajR, and thus the string tension 
aj is given by 



9^jU + 1) 



2Aj 



(1) 



In more a detailed fiux-tube model, the color-field Ej may vary in magnitude across the fiux- 
tube cross-section; the cross-section in Eq. (|l]) is then defined by Aj = (/ EjdAj)^/ J EjdAj] 
it is also possible to relax the constraint that Aj is independent of R. 

If the general features of the fiux-tube model of confinement are consistent with QCD, 
then the connection Eq. (|1]) between aj and Aj should hold. However, within the con- 
text of the fiux-tube picture, Aj is an unknown function of the quark representation. The 
cross-section is determined by the fundamental dynamics of QCD. Flux-tube formation has 
been observed in lattice QCD simulations in four dimensions 0J^,^,|,|],0 and as well as in 
three dimensions However, previous fiux-tube measurements have only been made for 
fundamental representation sources. 

This paper presents results of the first direct lattice measurements of Aj for static quarks 
in different representations of the gauge group. This work is aimed in part at establishing 
the connection Eq. (|l|) between the string tension and fiux-tube cross-section. As well, the 
results obtained here go beyond the flux-tube model, providing important new information 
about the dynamics underlying confinement. 

For example, if confinement is due to a bulk property of the QCD vacuum, such as a 
vacuum pressure, then Aj is expected to increase with representation (a natural "mechani- 
cal" response of a "medium" to the injection of more intense fields). Consequently, in such 
a scenario, the string tension aj is expected to increase with representation less rapidly than 
the quark Casimir. This situation is realized in a wide class of phenomenological models, 
including bag models ^j^, models based on a description of the QCD vacuum as a color- 



dielectric medium 10 , and some models of confinement based on vacuum condensation 
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The MIT bag ^ is typical of this general class of models. In the case of the heavy 
quark-antiquark {QjQj) system described above, the magnitude of the color-electric field is 
determined by a balance between the pressure generated by the field and an external "bag" 
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pressure B |]l|, ^Ej = B. It follows that Aj oc Qj, and the resulting string tension aj also 
scales as the square root of the Casimir, cxj oc (j(j + 1))^^^. 

In fact, the phenomenological models discussed above are not compatible with lattice 
simulations, which have shown that the string tension actually scales to an excellent ap- 
proximation like the Casimir of the representation. This has been observed for both SU(2) 
and SU(3) gauge groups ||T^. In the context of the flux-tube model, this suggests that the 



cross-section Aj is independent of representation. 

Direct measurements of the fiux-tube cross-section for heavy quark-antiquark [QjQj) 
sources are obtained here in three representations of three-dimensional SU(2) lattice gauge 
theory: fundamental (j = 1/2), adjoint (j = 1), and quartet (j = 3/2). A more thorough 
check of scaling and finite volume effects is achieved by working in three dimensions than 
would be obtained (with the same computing power) in four dimensions. We think that our 
results are relevant to the problem of confinement in four-dimensional QCD. In particular, 
previous lattice studies have shown that the string tension scales like the Casimir of the 



representation in three dimensions [^,0], as well as in four dimensions |T2[. Moreover, the 
fiux-tube picture of confinement is qualitatively the same in both three and four dimensions. 
It is therefore reasonable to expect that the qualitative features of flux-tubes reported here 
would be reproduced in four dimensions. Of course, such a calculation can and should be 
done. 

We find Aj ~ constant for the three representations, to within about 10% (a rough 
estimate of the overall quality of our data). This is consistent with the fiux-tube picture, 
given that the string tension scales like the Casimir of the representation (which is confirmed 
here up to the quartet representation). Several additional qualitative features of the fiux- 
tube picture are also verified. 

These results suggest a connection between confinement in QCD and the physics of 
a dual superconductor. Indeed, if a multiply- charged monople would be inserted into an 
ordinary (type II) superconductor, all the quanta of magnetic flux would be carried by a 
single flux-tube, whose diameter is fixed by the penetration depth A pair of monopoles 
of opposite sign would therefore be confined, with a string tension that would scale like the 
squared- charge. 

It is well-known that dual superconductivity (magnetic monopole condensation) results 
in confinement of electric charges in compact QED in three-dimensions (QED3) ||16|. A 
simple extension of the analytical calculation of Ref. [|T^] in the Villain approximation to 
the Wilson action, to include Wilson loops for multiply-charged sources, demonstrates that 
the string tension scales like the squared- charge. 

We have performed lattice simulations of singly- and doubly-charged Wilson loops in 
compact QED3, and our results confirm the expected scaling properties of the string tension 
and flux-tube cross-section. The potential is found to scale like the squared-charge to within 
a few percent, and the flux-tubes in the two cases have the same cross-section to within 
about 10%. The results of our three-dimensional SU(2) and U(l) simulations taken together 
lend some support, albeit indirectly, to the dual superconductor picture of confinement in 
four-dimensional QCD 
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II. METHOD 



To begin with, we consider the three-dimensional SU(2) lattice theory. Wilson loops 
are used to introduce static QjQj sources. Lattice measurements of the color-electric and 
-magnetic fields generated by these sources are obtained from correlators jFj"^ of plaquettes 
with a Wilson loop 



(2) 



where U^^{x) is the plaquette located at x (measured relative to the center of the Wilson 
loop), and Wj is the normalized trace of the Wilson loop in the j-th representation: 

W,^,^^A\[VAuA. (3) 



T>j\^i\ denotes an appropriate irreducible representation of the link f/;, and L the closed 
loop. [3 = 4:/{g'^a), where the coupling constant g has dimensions of (mass)^/^ in three 
dimensions. 

In the continuum limit, the trace of a 1 x 1 plaquette is by construction independent 
of representation (up to overall normalizations). As in several previous lattice calculations 
of higher representation Wilson loops (cf. Refs. 0,|l^,|l3), we use the action expressed in 
terms of links in the fundamental representation to perform simulations at arbitrary /3. The 
trace of the plaquette If^^ in the fundamental representation is also used to compute the 
correlators of Eq. (||) . 

In the continuum limit the correlator jFj"^ corresponds to the expectation value of the 
square of the Euchdean field strength F''^"' = d^'A'"' - S'^A"^ + ge^bcj^^bt^ji^cy ^ 

jta = - m^-nx (4) 

where the expectation value (. . ■)q .q, is taken in a state with external sources in the j-th 
representation, and (. . .)q is the vacuum expectation value. 

To compute the energy density, the Euclidean 3-axis is identified with a temporal side 
of the Wilson loop, and the 1-axis with a radial side. We separate contributions to the 
total energy density Stot corresponding to the two spatial components of the color-electric 
field (in the directions parallel and perpendicular to the line joining the quarks), and the 
color-magnetic field (a scalar in three dimensions): 

(x) = (x) + £^ (x) + 8f (x) , (5) 



where 



£\\{x) = -Tf{x), 

S^ix) ^ -^f (x), (6) 
Sf{x)^J^j\x). 
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Notice the relative minus sign between the electric and magnetic components of the Eu- 
clidean energy density. 

Previous flux measurements for fundamental representation sources have been made in 
four dimensions p|,^,^J^J^,0] , and in three dimensions f^. Following Haymaker and Wosiek 
, we achieve a significant enhancement in the signal to noise for the correlators by replacing 
Eq. (D with: 



(7) 



where xr is a reference point chosen far enough from the Wilson loop that the factorization 
{WjU{x{i)) ~ {Wj) {U) is satisfied. As in Ref. 0, we find that this happens well within 
the lattice volume. We verified explicitly that the right-hand side of Eq. is insensitive, 
within our statistical errors, to variations in xr over a wide range (when measurements are 
made for x in a region around the Wilson loop of sufficient size to suit our purposes). The 
results presented here were obtained with xr taken at a distance of half the lattice size from 
the center of the Wilson loop, in the direction transverse to the plane of the loop. The 
advantage to using Eq. (|^ is that the fluctuations in the product WjU^'^ , due mainly to 
the Wilson loop, tend to cancel in the vacuum subtraction when computed configuration by 
configuration. 

Another reduction in the statistical errors is readily achieved by performing some link 
integrations analytically, following the multihit procedure introduced by Parisi, Petronzio 
and Rapuano Consider a link variable Ui which appears linearly in the observable 



of interest. The simplest analytical integration over Ui takes account of nearest neighbor 
couplings in the action: 

Il[pKi) J 

where Ki is the sum of the four "staples" coupling to the link of interest Ui, and 

kiVi = Ki, detVi = l. (9) 



A further variance reduction has been developed by Mawhinney which takes account 
of effective next-to-nearest neighbor interactions with the link of interestlll 

where S is the action [only next-to-nearest neighbor couplings to Ui are relevant in Eq. (0)], 
and 



Mawhinney derived next-to-nearest neighbor variance reductions for fundamental and adjoint 
representations by employing an axial gauge-fixing [Q. We have generalized his result to arbitrary 
representations without gauge fixing. 



5 



KV/^j:Ui^{x)Y.Uif.i^ + l±)Ul{x), detV/^l. (11) 

i± fi^i± 

l± are four unit vectors perpendicular to /. The oriented plaquette ?7//i(x + l±) is computed 
with the hnk Ui{x + l±) appearing first on the left [x is the position of the base of the link 
Ui in Eq. (p!oD ]. The sum over unit vectors l± and ft in Eq. (|ll]) is taken over both parallel 



and antiparallel orientations with respect to a set of fixed basis vectors. An integration over 
the four links Ui{x + l±) is implicit in Eq. ([T0|). 

The second order variance reduction of Eq. ( |10|) cannot be applied to the links in the 
corners of the Wilson loop, since some links would then appear more than once in the 
integrand. Likewise, the first-order variance reduction Eq. (H) can only be applied to one 
link in a corner. Further restrictions apply to calculations of the plaquette correlators. 

For Wilson loops with less than six links on a side, we use the first-order variance 
reduction Eq. (^ for all links in the loop, except for one link at each corner, where no 
variance reduction is used. A plaquette correlator can be measured simultaneously provided 
that all sides of the plaquette are at least one node from the sides of the Wilson loop. 

For Wilson loops of size 6 x 6 or larger, we minimize the variance by using a combination 
of first- and second-order variance reductions. Equation (§) is appplied to the first link and 
to the second-to-last link on each (oriented) side of the Wilson loop; Eq. (0) is applied to 
all other links, except the last link on each side, where no variance reduction is used. In 
this case, a correlator can be measured simultaneously only if all sides of the plaquette are 
at least two nodes from the sides of the Wilson loop. [This variance reduction scheme can 
also be applied to Wilson loops as small as 4 x 4 if measurements of correlators near the 
center of the loop are not desired. This scheme is significantly more effective than the one 
employed in Ref. which uses only second-order variance reductions.] 



For large /3, /c/ ~ 4 and k'l ~ 12; Eq. (|T0]) then provides an estimate of the reduction v^ed 
in the variance of a Wilson loop using the above scheme, compared to the variance when 
only "unreduced" links are used (cf. Refs. pO|,p^ 



~ i h{m I \ /i(4,9) I ■ ^^^> 

For example, the variance reduction for a Wilson loop of size 6 x 6 in the quartet represen- 
tation at /? = 10 is estimated to be a factor of ~ 90. Our numerical results are consistent 
with Eq. (IID. 

The trace of an element of the group in the j-th representation can be expressed in terms 
of its trace in the fundamental representation using trigonometric relations among the group 



characters. In the case of the adjoint and quartet representations |21]: 



1^3/2 = 21^1% - 1^1/2, 

Wi = (mli^ - l) /3. (13) 

Hence one need only compute the Wilson loop in the fundamental representation, using the 
"unreduced" links or the "reduced" elements VJ, VJ' of Eqs. (^ and (pUj), as the case may 
be. The Wilson loops in higher representations then follow from Eq. ([T3|). The Bessel func- 
tions for the analytical integrations are tabulated separately for the three representations. 
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III. RESULTS AND DISCUSSION 



Our main results were obtained on a 32^ lattice at /5 = 10 (which is well within the 
scaling region for the string tension on a lattice of this size fl^). Wilson loops and plaquette 
correlators were calculated in the three representations j = 1/2, 1, and 3/2 for all loops of 
sizes T X R from 3 x 4 to 8 x 8 (these observables were measured in groups in several 
separate runs). Some additional data was taken at /? = 14 in order to check for scaling of 
the physical flux-tube dimensions. A standard heat-bath algorithm was employed. More 
than 10,000 sweeps were typically used for thermalization. 2,000 measurements were made, 
taking 20 sweeps between measurements. The resulting integrated autocorrelation times r^nt 
for the Wilson loops generally satisfy Tint ^ 1, consistent with the results of a systematic 
study made in Ref. [Q. Estimates of the statistical errors were obtained using the jackknife 
method. However, measurements of different observables (and of a given observable in the 
three representations) tend to be strongly correlated, since many Wilson loops and plaquette 
correlators were measured simultaneously on a given lattice. 

The quartet representation is much more difficult to measure than the two lower repre- 
sentations, due to the exponential suppression of the Wilson loop with the QjQj potential, 
which is found to scale with the Casimir of the representation. Energy density measurements 
in the quartet case obtained from loops larger than about 6x6 are of poor quality, although 
these data are consistent with conclusions drawn from results obtained from smaller loops. 

Representative data for Wilson loops in the three representations are shown in Fig. |1|. 
Earlier studies have shown that the potentials scale with the Casimir of the representation at 
essentially all lengths scales R [[l3ip!4| . This is made evident in Fig. |l|, where the logarithms 
of the Wilson loops are scaled by a ratio of Casimirs, 

' j(j + l) ^ ^ 

The quantity —\ia{Wj{T,R))/T, which extrapolates to the QjQj potential Vj{R) in the 
limit T —>■ oo, is found to scale as j{j + 1) to within a few tenths of a percent at all T 
and R considered here. A simple extrapolation of the data using Vj{R) ~ \n[{Wj{T^i,-^ — 
1, R)) / (W,(Tjnax5 where T^ax is the largest T value in the data set, gives agreement to 
a few tenths of a percent with the results of a careful statistical analysis of fundamental and 
adjoint Wilson loops reported in Ref. |]13 . 



Several attributes of the plaquette correlators were measured. To begin with, results for 
the fundamental and adjoint representations are presented. The correlators were measured 
over a range of distances x± from the center of the Wilson loop, in the direction normal 
to the plane of the loop. Results for the T x i? = 8 x 6 loop are shown in Fig. ^. The 
cross-sections of the fundamental and adjoint representation fiux-tubes are indistinguishable 
within statistical errors. This is true for all Wilson loops that were considered. For example, 
the T evolution of £']' for R = 6 Wilson loops is illustrated in Fig. ^. As observed in Refs. 
0J^, the plaquette correlators are more sensitive to higher states than the Wilson loop. Our 
data are consistent with a one-excited-state parameterization given in Ref. 

Figure ^ demonstrates that the component of the color-electric field parallel to the line 
joining the charges dominates the energy, as assumed in the fiux-tube model. The magnetic 
energy turns out to be negative, which has also been observed in four- dimensional SU(2) 
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lattice theory 0. The formation of a well-defined flux-tube is demonstrated by measure- 
ments of Sj in the plane of the Wilson loop. Figure H shows Sj for the T x R = 6 x 8 loop 
as a function of the longitudinal distance x\\ of the plaquette centroid from the center of 
the loop. Notice the approximate symmetry of the energy density about the center of the 
loop. The formation of the fiux-tube is further illustrated in Fig. where 8j is shown as a 
function of the radial separation R of the Wilson loop (for fixed T = 6) 

A stringent test of energy density calculations using Eq. is provided a sum rule 



derived by Michael p2 



a'E^r(^)=^.(^)- (15) 

The analogous sum rule in four- dimensional SU(2) was studied in detail by Haymaker and 
Woseik 1^. The flux-tube picture suggests a related sum rule that is much simpler to 
measure. If the interaction energy is dominated by a constant color-electric field along the 
line joining the charges (as expected in the limit of quark separations much greater than the 
fiux-tube thickness), then the integral of the energy density along one transverse "slice" of 
the fiux-tube should equal the string tension [cf. aj = lim^^oo(V,(-R) — Vj{R — a)) /a]: 

a ^ Sj{x±, X|| = fixed) ~ aj, (16) 

where the sum is taken over positive and negative distances x± from the plane of the Wilson 
loop. 

Our results are in good agreement with Eq. (|TB|). Figure |^ shows the left-hand-side of 
this equation for the T x R = 8 x 6 loop in the fundamental and adjoint representations, 
using a variable cutoff x*^ on the sum. The right-hand-side of Eq. (|TB]) is illustrated by the 
dashed line in Fig. M (our estimates of the string tension in the two representations agree 



with Ref. |Tj] to within a few percent). These results again demonstrate that the flux-tube 
cross-sections for the fundamental and adjoint repesentations are indistinguishable within 
statistical errors. 

Figures also demonstrate that the local energy densities scale to a good approxima- 
tion like the Casimir of the representation throughout the flux-tube. This is a very strong 
test of the validity of the flux-tube model for Aj = constant. Since the magnitude of the 
color-electric fleld varies across the flux-tube cross-section (cf. Fig. |^), a proper determi- 
nation of the numerical value of Aj should be made in terms of expectation values of the 
color-fleld, as described below Eq. (|1|) [some prescription for deflning the Abelian projection 
of the color-fleld would also be required]. However, a rough estimate Aj ~ 8a inferred from 
Fig. 1^ is consistent with Eq. (|l|), given the estimate of the string tension CjCPj ^ O.lig'^ (cf. 
Ref. H). 

The cross-section is also found to be approximately independent of R. A similar con- 
clusion was reached in four dimensions in Ref. H. On the other hand, the cross-section in 



the strong coupling limit in four dimensions is found to increase logarithmically with R |^ 
Within statistical errors the range in R considered here is not sufficient to rule out such a 
weak dependence on the radial separation. 

Our measurements of the quartet representation (j = 3/2) correlators are consistent with 
the above results. We compare data in the three representations taken from the T x R = 5x6 
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Wilson loop: energy density profiles transverse to the fiux-tube are shown in Fig. and 
the sum rules Eq. ([T6| ) in Fig. ^j. Data obtained from larger Wilson loops are consistent 
with these results although, as mentioned above, the quartet data for larger loops are of 
poor quahty, due to an exponential suppression of the Wilson loop with the Casimir of the 
representation. 

We checked for scaling of the physical fiux-tube dimensions by running at f3 = 14. The 
sum rule Eq. (|TB|) for the fundamental representation is compared at the two values of (3 in 
Fig. p. The cutoff x*j_ on the sum is expressed here in units of the physical coupling constant 
g. The data at /5 = 10 are for a 6 x 6 Wilson loop, while the data at /5 = 14 are for an 
8x8 loop. These Wilson loops have roughly the same dimensions in physical units (T and 
R measured in units of l/g"^)- 

These results show good evidence for scaling in the energy density and flux-tube cross- 
section (scaling is also observed in our adjoint and quartet representation data). However, 
the factorization assumed in Eq. (|^ breaks down in the /5 = 14 data at the largest values of 
the cutoff x]_ shown in Fig. ^ {x]_ ~ 9.5 in lattice units, to be compared with xr = 16); the 
sum is found to diverge linearly with x*j_ at large cutoffs. A similar behavior was observed 
in four-dimensional lattice calculations in Ref. 0, where a correction for this effect was 
proposed. Nevertheless, scaling of the cross-section is clearly supported by data in the 
region x]_ < 1-8/5^^. 

As described in the Introduction, the results of our SU(2) simulation suggest a connection 
between conflnement in QCD and the physics of a dual superconductor. In this connection, 
we have calculated Wilson loops Wn in compact QED3 for singly- and doubly-charged sources 
(n = l,2): 

W„ = Re liiUiT, (17) 

where the phase Ui for the link I deflnes the singly-charged representation (i.e., Ui is the 
phase used to compute the Wilson action). The string tension is expected to scale like 
the squared-charge, as demonstrated by an extension of the Villain approximation used in 
Ref. [jl^ to include multiply-charged Wilson loops. As in an ordinary superconductor, the 



flux-tube cross-section is expected to be independent of source charge. 

Wilson loops and plaquette correlators were measured on a 32^ lattice at /3 = 2.4. More 
than 10,000 sweeps were used for thermalization, and 1,000 measurements were made (90 
sweeps were taken between measurements). Variance reduction methods similar to those 
used in our SU(2) simulations were employed.H 

Results for the Wilson loops are given in Fig. |1^. Some data for the triply-charged 
Wilson loop {n = 3) are also shown (useful measurements of the plaquette correlators for 
n = 3 would require much larger statistics). Estimates of the potential for n = 1 obtained 
from a simple extrapolation of these data are consistent with results presented in graphical 
form in Ref. |E3|1. 



^The analytical integrals given in Eqs. (^) and (10) are easily adapted to the U(l) theory. 
The SU(2) link VjlUi] becomes {Ui)"^, the Bessel function ratios l2j+i{x)/ Ii{x) are replaced by 



In{x)/Io{x) 



and det(y) becomes ahs{V). 
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We find that — In {Wn{R,T)) /T, which extrapolates to the potential Vn{R) in the limit 
T — > cxo, scales like to within about 2% for all T and R that were considered, in good 
agreement with the expected scaling properties of the string tension. However, the deviation 
from scaling is about an order of magnitude larger than the statistical errors in the data. 
String vibrational modes are known to make a significant contribution to Vi in the range 
of R considered here (lattice simulations pSf are in agreement with theoretical expectations 
2^). Simple arguments |]T3[ suggest that the vibrational term in Vn may scale like n, which 



could account for the small deviation from n"^ scaling in the logarithms of the Wilson loops. 

The energy sum rule analogous to Eq. (ITBp for the T x R = 5 x 5 loop is shown in Fig. 
|TT]. The dashed line shows the n = 1 string tension taken from Ref. ||2^. These results 
provide the first direct evidence from lattice simulations that the fiux-tube cross-section in 
compact QED3 is independent of the source charge, the expected behavior in the case of a 
(dual) superconducting medium. 



IV. SUMMARY 

The first direct measurements of the fiux-tube cross-section as a function of representa- 
tion in SU(2) lattice gauge theory were made. We found Aj ~ constant, to within about 10% 
(a rough estimate of the overall quality of our data) for the three representations j = 1/2, 
1, and 3/2. Our results are consistent with a connection between the string tension and 
cross-section suggested by a simplified model of fiux-tube formation, aj = g'^j{j + l)/{2Aj), 
given that the string tension scales like the Casimir j{j + 1), as observed in previous lattice 
studies in both three and four dimensions (and confirmed here up to the quartet represen- 
tation). We also confirmed several additional quahtative features of the fiux-tube picture of 
color-electric confinement. These results can be used discriminate among phenomenological 
models of the physics underlying confinement. For example, many models in which confine- 
ment is due to a bulk property of the QCD vacuum (such as a vacuum pressure) predict a 
sufficiently rapid increase in Aj with representation as to be incompatible with the results 
obtained from our lattice simulations. 

We also made fiux-tube measurements in compact QED3, which exhibits electric confine- 
ment due to magnetic monopole condensation. We considered singly- and doubly-charged 
Wilson loops. The string tension was found to scale like the squared- charge, and the fiux- 
tube cross-section was found to be independent of the charge, to a good approximation. 
The results of our three-dimensional SU(2) and U(l) simulations taken together lend some 
support, albeit indirectly, to a conjecture that the dual superconductor mechanism underlies 
confinement in compact gauge theories in both three and four dimensions. This conclusion 
is also supported by the results of a recent study of dual Abrikosov vortices in an Abelian 
projection of SU(2) lattice gauge theory in four dimensions |2^. Flux-tube measurements 



in four-dimensional SU(2) gauge theory similar to those reported here should be made in 
order to further explore this possibility. 
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FIGURES 



FIG. 1. T evolution of Wilson loops in three representations of SU(2) lattice gauge theory: 
j = 1/2 (o), j = 1 (□), and j = 3/2 (A). Cj is a ratio of Casimirs, defined in Eq. (p^). The 
quantity —\n(Wj{T^R))/T extrapolates to the QjQj potential in the limit T — > oo. 

FIG. 2. Energy density profiles transverse to the plane of the T x R = 8x6 Wilson loop [j = 1/2 
(o), and j = 1 (□)]. These results are an average over plaquettes with centroids at distances ±x± 
transverse to the plane of the loop. 

FIG. 3. T evolution of = 0, 2a) for Wilson loops with R = 6 [j = 1/2 (o), j = 1 (□)]. 

FIG. 4. Energy density profile in the plane of the T x i? = 6 x 8 Wilson loop [j = 1/2 (o), 
j = 1 (n)]. a^ii is the distance of the centroid of the plaquette from the center of the Wilson loop. 
The radial sides of the Wilson loop are located at xy = ib4a. Plaquettes with a side touching the 
Wilson loop cannot be measured using the variance reduction of Eq. and are not shown. 

FIG. 5. Energy density £j{x± = 0, 2a) as a function of R, for fixed T = 6 [j = 1/2 (o), j = 1 
(□)]■ 



FIG. 6. Energy sum rule Eq. (|T|) for the T x = 8 x 6 Wilson loop [j = 1/2 (o), j = 1 (□)]. 
The sum in Eq. (16) is evaluated using a variable cutoff x*j_. The dashed line shows the scaled 
string tension cjaj, estimated to about 5% (cf. Ref. 0]). 



FIG. 7. Energy density profile transverse to the plane of the T x R = 5 x 6 Wilson loop, in the 
fundamental and quartet representations \j = 1/2 (o), j = 3/2 (A)]. 



FIG. 8. Energy sum rule Eq. (|lq ) for the T xR = 5x6 Wilson loop, in all three representations 
[j = 1/2 (o), j = 1 (□), j = 3/2 (A)]. The sum rule improves with increasing T and R (cf. Fig. 

!)• 



FIG. 9. Scaling of the energy sum rule Eq. (16) for the fundamental representation. The open 
data points were taken at /3 = 10 (for a 6 x 6 Wilson loop), and the filled points at /3 = 14 (for an 
8x8 loop). The cutoff x*j_ is expressed here in units of the physical coupling constant g. 



FIG. 10. T evolution of multiply-charged Wilson loops in compact QED3 [n = 1 (o), n = 2 
(□), and n = 3 (A)]. 
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FIG. 11. Energy sum rule Eq. ( |l6|) for singly- and doubly-charged sources in compact QED3, 
for the T X i? = 5 X 5 Wilson loop [n = 1 (o), n = 2 (□)]. The dashed line shows the string tension 
cr„=i, estimated to about 10% in Ref. p3| . 
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